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We study the formation and final structure of the rare earth peak (A ~ 160) of the r-process 
nucleosynthesis. The rare earth peak forms at late times in the r-process after neutron exhaustion 
(neutron-to-seed ratio unity or R — 1) as matter decays back to stability. Since rare earth peak 
formation does not occur during (n, 7) (7, n) equilibrium it is sensitive to the strong interplay 
between late time thermodynamic evolution and nuclear physics input. Depending on the conditions 
the peak forms either because of the pattern of the neutron capture rates or because of the pattern of 
the separation energies. We analyze three mass models under different thermodynamic conditions. 
We find that the subtleties of each mass model, including separation energies and neutron capture 
rates, influence not only the final shape of the peak but also when it forms. We identify the range 
of nuclei which are influential in rare earth peak formation. 

PACS numbers: 20.30.-k, 26.30.Hj, 26.50.+X 



I. INTRODUCTION 

Approximately half of the elements beyond A > 100 
are made in the 'rapid' neutron capture process, or re- 
process, in which successive neutron captures occur on 
timescales faster than /3-decays. At the present time, 
there is significant uncertainty with the astrophysical en- 
vironment responsible for this synthesis event [I] [2] . The 
leading candidate site [3] is believed to be core-collapse 
supernovae e.g [HUP] even though most recent simula- 
tions do not yield favorable conditions for the r-process 
[T1TH3] . Other candidate sites include compact object 
mergers [TH - f2"T] , gamma-ray burst outflows [2"2Tf2~i] , neu- 
trino induced nucleosynthesis in He shells [351 12S] , super- 
nova fallback [57], and collapse of O-Ne-Mg cores [2"5rl50"] . 

Experimentally, it is difficult to measure the properties 
of the short-lived nuclei far from stability that participate 
in the r-process. Recent developments using radioactive 
beams show promise (e.g. |31l 132] ). but current exper- 
imental data on neutron-rich isotopes is limited. Thus 
r-process studies must rely not only on model calcula- 
tions of the environment, but also on theoretical mass 
models, e.g. [53TI55] . 

Despite these difficulties, much has been learned about 
the r-process over the past 50 years. The most prominent 
features in the r-process abundance distribution above 
atomic mass number of A = 100 are two distinct peaks 
occurring at A = 130 and A = 195. It was hypothesized 
very early that the formation of these peaks should be 
associated with the long /3-decay rates of closed neutron 
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shells p86l . Since this seminal paper much effort has been 
put in to researching the conditions for a sufficient initial 
neutron-to-seed ratio, a key requirement in order to pro- 
duce a 'main' r-process out to the third peak (A = 195). 
For reviews see [T] 12] 157] 155] . 

After the two main peaks, the second most prominent 
feature above A=100 is the smaller peak near A ~ 160 
known as the rare earth peak. While less abundant than 
the other peaks, the rare earth peak can in principle be 
used as a powerful tool and offers an alternative way to 
probe the r-process. This is due to the following proper- 
ties: (1) Observational data from metal-poor stars show 
very consistent trends among the rare earth and heavier 
elements. This suggests that these elements were cre- 
ated in the same type of synthesis event [55]. Thus, 
the rare earth peak provides a natural diagnostic of re- 
process models. (2) The rare earth peak forms away from 
closed neutron or proton shells and therefore by a differ- 
ent mechanism than the main peaks. This means it is 
a different and unique probe of late-time r-process con- 
ditions. (3) The rare earth peak is extremely sensitive 
not only to late-time thermodynamic behavior, but also 
to nuclear physics input [iOTfl2"] . Typical variations in fi- 
nal rare earth abundance patterns from simulations with 
different mass models are highlighted in Figure [T] 

To date the rare earth region has received relatively 
little attention. Fission cycling has been suggested as 
a mechanism for obtaining the rare earth peak [331 H3] , 
but it is not favored [45] . Large uncertainties found in 
fission probabilities and fragment distributions of current 
nuclear models further compound difficulties with a suc- 
cessful description of rare earth peak formation by fission 
cycling JT5J- Surman et al. [ID] investigated the forma- 
tion of the peak in a hot r-process environment with tem- 
peratures high enough to support (n, 7) ^ (7, n) equilib- 
rium. The formation of the rare earth peak under these 
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conditions was attributed to the co-action of nuclear de- 
formation and /3-decay as the free neutrons are quickly 
captured during freeze-out. This was followed by a study 
of late-time abundances changes among the major peaks 
[47S] , Otsuki et al. [37J investigated a range of r-process 
models and found similar rare earth elemental abundance 
patterns, provided the temperature was constant during 
freeze-out. Most recently, Arcones et al. [31] studied the 
sensitivity of late-time abundance fluctuations to changes 
in the nuclear physics inputs. Arcones et al. [42] pointed 
out that the rare earth peak is sensitive to changes at 
late-times, e.g. to non-equilibrium effects such as neu- 
tron capture even when the abundance of free neutrons 
can become very low (~ 10~ 5 ). 

This manuscript presents a more complete picture of 
rare earth peak formation. We explore the sensitivity 
of the peak formation mechanism to late-time thermo- 
dynamic behavior and nuclear physics input. The 'fun- 
neling' formation mechanism of [40] is reviewed for hot 
evolutions. We introduce a different 'trapping' mech- 
anism for peak formation in cold evolutions where the 
temperatures and densities decline relatively quickly and 
therefore photo-dissociation plays no role in the late-time 
dynamics after R = 1. We study the effects of three dif- 
ferent mass models and show how large uncertainties in 
this region stem from nuclear physics. Lastly, we show 
that the nuclei which contribute to peak formation are 
approximately 10 to 15 neutrons from stability, and thus 
represent prime candidates to be measured in future ra- 
dioactive ion beam facilities (FRIB gS] or FAIR [33 )■ 



II. R-PROCESS CONDITIONS AND 
CALCULATIONS 

Abundance weighted lifetimes are used throughout the 
text to characterize the late-time dynamics of the re- 
process. These are provided below for the reader's con- 
venience: 
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where N n is the neutron number density, (av)z.A the 
thermally averaged neutron capture cross section for nu- 
clei (Z,A), X in (Z,A) the photo-dissociation rate for nu- 
clei {Z,A), \{i(Z,A) the full /3-decay rate (including /3- 
delayed neutron emission channels) for nuclei {Z 1 A) and 
Y(Z,A) the abundance of nuclei (Z,A). A reduced sum 



denoted with a superscript "REP" is taken over the rare 
earth region, A = 150 to A = 180, when applicable. The 
neutron-to-seed ratio or R is defined as: 
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where Y n is the abundance of free neutrons. 

Since rare earth peak formation is highly dependent 
on the rate of decrease in the temperature and density, 
we consider rare earth peak formation under two differ- 
ent thermodynamic evolutions. One scenario is a clas- 
sical 'hot' r-process which operates under high temper- 
atures (Tg > 1) at the time in which neutron captures 
are important for peak formation. A second scenario is a 
'cold' r-process which operates under low temperatures 
(Tg ~ 0.5) at the time in which neutron captures are 
important for peak formation |50j . 

The classical r-process begins with a phase of (n, 7) ^ 
(7, n) equilibrium marked by an abundance weighted life- 
time ratio of neutron capture to photo-dissociation of 
t„ 7 /t 7JI = 1. During this phase the temperature is still 
sufficiently high so that neutron captures dominate /3- 
decays (t^/t,^^ 1) and the Saha equation can be used 
to determine abundances along an isotopic chain [55] , 

The second phase, known as the freeze-out epoch, is 
marked by the weakening of the (n,7) ^ equi- 
librium (T„ 7 /T 7Ji < 1) and the abundance weighted life- 
time ratio of /3-decay versus neutron capture falls to 
Tp/r n y~ 1. It is during this phase that the formation of 
the rare earth peak proceeds with competition between 
neutron captures, photo-disintegrations and /3-decays. 

In the cold r-process the first phase (71,7) <=± {j,n) 
equilibrium is dramatically shorter than the first phase of 
the classical scenario. Freeze-out is now caused by a rapid 
drop in temperature rather than the consumption of free 
neutrons (as in the classical case). The bulk of the cold 
r-process operates in the second phase, under low tem- 
peratures (Tg ~ 0.5), where photo-disintegrations have 
frozen out [50] . 

Once neutron exhaustion (R ~ 1) occurs in the cold 
r-process the free neutrons available to the system must 
come from the recapture of /3-delayed emitted neutrons. 
The importance of this effect on the final abundance 
distribution was noted in [42] [51]. This recapture ef- 
fect is crucial to peak formation as can be seen from 
the fact that malformed abundance distributions result 
if /3-delayed neutron emission is artificially turned off (see 
[42]). 

Our calculations consists of a nuclear reaction net- 
work containing r-process relevant nuclides as described 
in [301 EE] ■ Previous versions of this network code have 
been used in the studies of Beun [S5] and Surman [55] , 
The primary reaction channels for nuclides in this section 
of the reaction network are beta-decay, neutron capture, 
and photo-dissociation. Our fully implicit r-process reac- 
tion network handles consistently neutron capture rates 
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at low temperatures and calculations with low abun- 
dances of free neutrons, both important for simulations 
with cold evolutions. For the initial abundances we use 
self-consistent output from an intermediate reaction net- 
work [SI] with PARDISO solver [55]. 

Our r-process calculations start at Tg = 2 with densi- 
ties p w 9 • 10 s g/cm 3 for hot evolutions with Y e = .30 
and p « 5 • 10 8 g/cm 3 for cold evolutions with Y e = .40. 
At this time the neutron-to-seed ratios are R « 45 and 
R w 35 respectively. We study the late-time hot and 
cold r-process evolutions in the context of a monotoni- 
cally decreasing temperature with density parameterized 
as: 

p(t) oc t~ n (3) 

where n controls the type of late-time r-process evolution 
(the time when rare earth peak formation occurs). For 
hot r-process evolutions we set n = 2 and for cold re- 
process evolutions we set n = 6. A decaying density of 
n = 2 is characteristic of wind models [56] [57] at late 
times while n — 6 represents a faster decline. 

We use three different mass models in our nucleosyn- 
thesis calculations: Finite Range Droplet Model (FRDM) 
[33j . Extended Thomas- Fermi with Strutinsky Integral 
and Quenching (ETFSI-Q) (3J and version 17 of the 
Hartree Fock Bogoliubov masses (HFB-17) 35 . The 
FRDM and ETFSI-Q neutron capture rates are from 
[58j and were computed with the statistical model code 
NON-SMOKER [5J]. The HFB-17 neutron capture rates 
are from the publicly available Brusslib online-database 
[60j and were computed with the statistical model code 
TALYS [3], which is also publicly available. The HFB 
mass model is under constant development and is there- 
fore updated with the latest experimental data and the- 
oretical techniques [52]. The /3-decay rates used in our 
r-process network come from [63]. 

III. PEAK FORMATION IN HOT 
ENVIRONMENTS 

The mechanism for rare earth peak formation in hot 
environments was first described in [ID]. We review the 
basic physical arguments in this section. 

Under hot conditions the r-process path (time ordered 
set of most abundant isotopes) traverses the NZ-plane 
between the line of stability and the neutron drip line. 
The path is initially constrained by (71,7) <=^ (7, ri) equi- 
librium and is thus found to lie on a line of constant sepa- 
ration energy via the Saha equation. As the free neutrons 
are consumed, the path moves back toward stability and 
(n, 7) <^ (7, n) equilibrium begins to break down. During 
this freeze-out from equilibrium, rare earth peak forma- 
tion can potentially occur. 

The necessary and sufficient conditions for peak forma- 
tion are as follows: (1) a deformation maximum or other 
nuclear structure effect must produce a kink in the lines 
of constant neutron separation energy around A ~ 160, 



and (2) the r-process path must traverse this kink re- 
gion during freeze-out, before /3-decay takes over in the 
region. The latter allows for the interplay of neutron 
capture, photo-dissociation and /3-decay as the r-process 
path crosses the region which contains the separation en- 
ergy kink. 

During peak formation, the r-process path moves to- 
ward stability at a rate approximately equal to the aver- 
age /3-decay rate along the path. The separation energy 
kink causes a corresponding kink in the r-process path 
as material moves through this important region. This 
provides a mismatch between the /3-decay rates of mate- 
rial below and above the kink. Due to the kink in the 
path, nuclei below the peak (A = 150 to A = 158) are 
farther from stability and so /3-decay faster than the av- 
erage nuclei along the path. Since the nuclei below the 
peak decay faster than the path moves, these nuclei then 
proceed to capture neutrons in an attempt to return the 
r-process path back to equilibrium. Conversely, due to 
the kink in the path, nuclei above the peak (A — 168 to 
A = 180) are closer to stability and so /3-decay slower 
than average along the path. The path therefore moves 
before these nuclei have a chance to decay and so they 
photo-dissociate to shift the r-process path back to equi- 
librium. In the peak region (A = 159 to A = 167) some 
nuclei are still in (n, 7) <=2 (7, n) equilibrium which limits 
the amount of material flowing out of the peak region in 
either direction. The net result causes material to funnel 
into the peak region, creating the local maximum. 

The essence of this effect is shown in Figure [2] At 
neutron exhaustion, R = 1 (left panel), the r-process is 
just beginning to break from (n, 7) <=^ (7, n) equilibrium. 
Here the path lies along a line of constant separation 
energy (~ 3.0 MeV) and the abundances show an odd- 
even effect due to the population of primarily even-N 
nuclei in equilibrium. No peak exists at this time. 

Later in the simulation (right panel), peak formation 
occurs as the path encounters the region with the sep- 
aration energy kink. The separation energy kink causes 
the kink in the r-process path. Nuclei along the path in 
the peak region have /3-decay rates which range from 1 
s _1 (above the kink) to 10 s _1 (below the kink). The 
resultant photo-dissociation above the kink and /3-decay 
followed by neutron capture below the kink causes mate- 
rial to funnel into the peak region. 



IV. PEAK FORMATION IN COLD 
ENVIRONMENTS 

In the previous section we analyzed rare earth peak 
formation in hot evolutions and found that photo- 
dissociation was crucial in peak formation. However, we 
also find well formed solar-like rare earth peaks in sim- 
ulations of cold environments where photo-dissociation 
plays no role in the dynamics after R = 1. 

After R = 1, the cold r-process path is controlled on 
average by the competition between neutron captures 
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and /3-decays 7-/3 /t„ 7 « 1. Locally, over the rare earth re- 
gion, the exact position of the path is more complicated 
due to the variation among individual rates. 

As the material decays back to stability peak formation 
will ensue if the path encounters a peak region where neu- 
tron capture rates are slow relative to the above and be- 
low regions. The essence of the effect is that slow neutron 
capture rates in the peak region cause a bow (inwards to- 
wards stability) in the lines of constant neutron capture 
rates relative to the lines of constant /3-decay rates thus 
causing material to become trapped in the peak region. 

The cold formation mechanism is shown in Figure [3] 
The left panel shows a snapshot of the abundance pattern 
and rates at neutron exhaustion, R = 1. At this point 
in time the r-process path is still influenced by residual 
photo-dissociation flows. This is reflected in an odd-even 
effect in the abundances and flat r-process path (simi- 
lar to hot evolutions). However, the photo-dissociation 
rates are decreasing so rapidly they play no further role in 
the dynamics after this point. Shortly, the neutron cap- 
ture rates will become comparable to the /3-decay rates 
and large odd-even behavior of the abundances will be 
washed-out [53] ■ In fact, this has already begun to hap- 
pen as can be seen with the slight bowing of the neu- 
tron capture rate lines in the peak region (A = 159 to 
A = 167). 

At a slightly later time in the simulation (right panel of 
Figure[3]) the system has moved closer to stability and the 
r-process path now encounters the slower capture rates 
in the peak region. Below the peak (A = 150 to A = 158) 
neutron captures occur much faster than /3-decay rates 
along the r-process path, so the net result is material 
shifting towards the peak region. In the peak region the 
path encounters the slow capture rates (note the bowing 
of the neutron capture rate lines) so that any material be- 
ing shifted into the peak region becomes hung up. Above 
the peak (A — 168 to A = 180) the flow of material is 
again dominated by the relatively faster neutron capture 
rates. The net result is trapping of material into the peak 
region. 

Another interesting feature in the right panel of Figure 
[3] is the trough to the left of the peak. A trough can 
occur if a gap in the r-process path proceeds for long 
periods of time as matter decays back to stability. Along 
a gap in the r-process path the neutron capture rates are 
relatively fast resulting in movement of material to more 
neutron-rich isotopes and a depletion of material in the 
gap region. 

In our figures, the lines of constant neutron capture 
rates have been averaged over even-N nuclei. Even-N 
neutron capture rates are more important to rare earth 
peak formation because at a given temperature, odd-N 
nuclei have faster neutron capture rates which causes ma- 
terial to pass through the odd-N nuclei quickly. Thus 
material builds up (or stays) in even-N nuclei which sets 
the r-process path. The importance of individual neu- 
tron capture rates in the rare earth peak was highlighted 
in [H]. 



V. INFLUENCE OF MASS MODEL ON RARE 
EARTH PEAK FORMATION 

From the previous two sections it is clear that the de- 
tails of the late-time thermodynamic evolution are crit- 
ical in setting the relevant nuclear physics and thus de- 
termine the mechanism for peak formation. 

Despite the differences in peak formation mechanisms, 
we find that the final abundances among simulations with 
the same mass model yet differing late-time thermody- 
namic behavior can be remarkably similar. This is in 
contrast to the differences found in the final abundance 
pattern when comparing between mass models with sim- 
ilar thermodynamic conditions. In this section we focus 
on the influence of different separation energies and neu- 
tron capture rates on rare earth peak formation. 

A successful peak formation is imprinted on the fi- 
nal abundances in a cold evolution when the r-process 
path encounters structure in the neutron capture rates 
and this structure lasts until the point at which /3- 
decays take over neutron captures in the region (t pep ss 
a few t pep ). 

A successful peak formation occurs in a hot evolution 
when the r-process path encounters a local deformation 
maximum leading to a well-defined kink structure in the 
separation energies in the rare earth region. 

For a given mass model, the structure of neutron cap- 
ture rates and the structure of the separation energies 
may not align in the NZ-plane. This in turn can affect 
the timing and location of peak formation and hence the 
nuclei which are relevant. 

Odd-even effects in the abundances can accumulate or 
persist through the decay back to stability resulting in 
visible features in the final abundances. Smoothing of 
the abundances typically occurs in between neutron cap- 
ture freeze-out (t pep = rj* ) and the time in which 
/3-decays fully take over neutron captures in the region 
(r£f p « afewrf p ). 

We now discuss three different mass models in this con- 
text. Since separation energies vary among mass models 
we instead (for consistency) use (SN), the abundance 
weighted average neutrons from stability, to measure the 
r-process path's progression. 

Compared to the other mass models studied here, we 
find that simulations which use the FRDM mass model 
best match the solar data in the rare earth peak region in 
both hot and cold evolutions. In fact we find (in agree- 
ment with previous studies [301112]) that the FRDM mass 
model is the only model to show a well-formed rare earth 
peak consistently in the final abundance pattern. 

Simulations with the FRDM mass model do not consis- 
tently form rare earth peaks far from stability ((SN) > 
20). Instead, peak formation ensues when the path is 
much closer; on average in between 15 and 20 neutrons 
away from stability. We can see the evolution of the 
peak region for a cold FRDM evolution in Figure [3] At 
(SN) ~ 20 (top panel) the structure in the capture rates 
has yet to manifest itself resulting in relatively flat abun- 
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dances. As the path moves back to stability, (SN) ~ 15 
(middle panel), it encounters nuclei in the peak region 
with relatively slower neutron capture rates than the 
surrounding regions (note the bending in the red lines). 
These conditions persist all the way back to stability re- 
sulting in a well-formed rare earth peak. A similar sce- 
nario occurs in hot evolutions; see Figure [5] 

FRDM shows a slight overlap between neutron cap- 
ture structure and separation energy structure. The 
structure in the separation energies occurring between 
(SN) ~ 12 to 20 and the structure in the capture rates 
occurring between (SN) ~ 10 to 15. This delays peak 
formation in cold scenarios until around 15 neutrons from 
stability, while hot evolutions typically begin peak forma- 
tion approximately 20 neutrons from stability. 

Simulations with the ETFSI mass model consistently 
form a solar-like rare earth peak far from the stable nuclei 
((SN) > 20). This is most apparent in colder simulations 
(see right panel of Figure|3]). However, this is not the end 
of the story as the material must decay back to stabil- 
ity. Figure [6] highlights this transition at an abundance 
weighted average of (SN) ~ 20 (top panel), 15 (middle 
panel) and 10 (bottom panel) neutrons from stability. As 
the decay back to stability proceeds the r-process path 
encounters nuclei whose neutron capture rates become 
homogeneous around the peak region. This slowly dis- 
solves the structure, flattening the lines of constant neu- 
tron capture rates (compare top and middle panels). By 
the time the path is on average 15 neutrons away from 
stability (middle panel) the cold trapping mechanism can 
not continue because neutron capture rates in the peak 
region are no longer slower than the surrounding regions. 
These conditions persist back to stability resulting in a 
final abundance pattern with a more modest rare earth 
peak. Note that a small odd-even effect reappears since 
in this model /3-delayed neutron emission is still relevant 
as neutron capture freezes out in the rare earth region 

(t-REP _ T REP\ 

Solar-like rare earth peaks form far from the stable nu- 
clei in ETFSI models under hot evolutions as well. Far 
from stability, the structure (kink) in the separation en- 
ergies results in the hot peak formation mechanism. Like 
the FRDM case, the separation energy kink in ETFSI 
disappears as one moves closer to stability. However, the 
kink disappears while neutron captures are still dominant 
( T nj P £ T p EP ) far lrom stability ({SN) > 20) resulting 
in a flattened final abundance distribution; see Figure [7] 

In this mass model the structure in the separation en- 
ergies occurs farther from stability {(SN) ^ 20) than the 
structure seen in the neutron capture rates ((SN) ~ 20) 
influencing peak formation in a similar fashion to the 
FRDM case. The gross separation energy structure oc- 
curs very early on "before" the top panel of Figure [7] and 
has already dissolved by (SN) ~ 20. 

Version 17 of the HFB mass model is optimized to 
over 2000 measured masses from |65j corresponding to 
a root mean square error of < 0.6MeV. This data set 
features detailed structure in the separation energies but 



little overall structure in the neutron capture rates for 
the nuclei relevant to rare earth peak formation. These 
features are reflected in our r-process abundances. 

Figure [8] shows the decay back to stability of a cold 
r-process using HFB-17. At every snapshot, highlight- 
ing the r-process path's decay back to stability, we do 
not find the structure in the neutron capture rates as is 
found in the other two mass models. It is this relative 
homogeneity in the neutron capture rates throughout the 
rare earth region which prevents the trapping mechanism 
from occurring in cold evolutions. 

In hot r-process evolutions the situation is more in- 
tricate than for the corresponding cases of the other two 
mass models. The detailed structure in the separation en- 
ergies results in a complex separation energy kink struc- 
ture in the rare earth region. However, due to the lack 
of gross structure as the separation energy increases (i.e. 
during the decay back to stability) the funncling mecha- 
nism cannot operate. This can be seen in Figure [9] and il- 
lustrates the subtleties involved in forming the rare earth 
peak. 

The discussion in this section showcases the need for 
nuclear structure measurements far from stability. As we 
have seen, the nuclei that are important for rare earth 
peak formation lie in between 10 and 20 neutrons away 
from stability. Furthermore, it is the nuclei which are the 
closest to stability, those in between 10 and 15 neutrons 
from stability, which are most influential to peak forma- 
tion as they set or potentially dissolve the peak structure 
all together. In Figure [10] we highlight these influential 
nuclei together with recent experimental mass measure- 
ments ( |66| [6 7) green and |65j gray) and known neutron 
capture rates ([53] red). 



VI. SUMMARY AND CONCLUSIONS 

We have studied the formation and evolution of the 
rare earth peak at late-times during the r-process. To 
take into account uncertainties with nuclear physics in 
the region our calculations employed three mass models 
(FRDM, ETFSI, and HFB-17). 

Two late-time evolutions were considered: A hot 
r-process with temperatures high enough to support 
(n, 7) (7, n) equilibrium and a cold r-process with 
lower temperatures where there are no photo-dissociation 
flows, only competition between neutron captures and j3- 
decays after R = 1. Both of these evolutions are similar 
at early times so that the changes in abundances are not 
due to to the physics that sets the neutron-to-seed ra- 
tio, but instead due to the changes in the nuclear physics 
input or changes in the late-time behavior of the evolu- 
tion. The differences in late-time evolution (hot vs cold) 
determine which nuclear physics input is important (sep- 
aration energies vs neutron capture rates respectively) 
during the final stages of the r-process. 

In hot evolutions the combination of photo- 
dissociation, neutron capture and beta-decay results 
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in a mechanism which funnels material into the peak 
region. A successful peak formation in hot evolutions is 
imprinted on the abundance pattern when the structure 
in the separation energies, the 'kink', is well defined and 
the r-process path crosses the kink region during the 
(n, 7) <^ (7,Ji) freeze-out. 

We contrast this with the peak formation mechanism 
which occurs in cold r-process environments. Here the 
important nuclear physics for peak formation lies in the 
local structure of the neutron capture rates. When the 
neutron capture rates are slow in the peak region rela- 
tive to the surrounding regions (creating the characteris- 
tic 'bow' in the lines of constant neutron capture rates) 
material can become trapped in the peak region, thus 
forming the peak. A successful peak formation in cold 
evolutions is imprinted on the abundance pattern when 
the structure in the neutron capture rates lasts until the 
point at which /3-decays take over neutron captures in 
the region (t^ p w a few t pep ). 

The rare earth peak is extremely sensitive to the sub- 
tleties of nuclear physics input. Neutron capture is par- 
ticularly important in both hot and cold evolutions. For 
instance, we find that neutron capture can play two com- 
peting roles in peak formation: it can be responsible for 
creating the peak, but also for potentially dissolving the 
peak (wash-out). Neutron capture rate structure and 



separation energy structure in the same mass model may 
not overlap in the NZ-plane. This in turn can affect the 
timing and location of peak formation in different ther- 
modynamic conditions. 

We have shown that the rare earth peak in principle 
offers unique insight into the late-time behavior of the 
r-process because it forms away from the closed shells 
during freeze-out while material decays back to stability. 
Rare earth peak formation is sensitive to the structure of 
separation energies and / or neutron capture rates about 
10 to 15 neutrons away from the stable rare earth peak. 
Future measurements at radioactive ion beam facilities 
should reach this important region and will be critical in 
placing constraints on nuclear models. This in turn will 
lead to improved r-process predictions; allowing the rare 
earth peak to evolve into a powerful tool for understand- 
ing the r-process. 
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FIG. 1: The resultant rare earth peaks from simulations with different mass models. FRDM (triangles), ETFSI 

[33] (squares) and HFB-17 [35] (circles) are shown along with the solar r-process abundance pattern, black line, 
Nq jT versus atomic mass (data from |69j). The same colors and geometric markers for each mass model will be used 

in the remaining figures. 
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FIG. 2: Shows how the rare earth peak forms under hot environments when the r-process path encounters the kink 

in the separation energies. Each left and right panel shows a snapshot of an abundance pattern along with 
separation energies (MeV), /3-decay rates (s -1 ), and r-process path from a simulation with the FRDM mass model. 

The left panel is a snapshot at neutron exhaustion, R = 1. The right panel is a snapshot of the r-process path 
beginning to move closer to stability while the peak forms. Connected black triangles represent the r-process path. 
Blue lines represent constant total /3-decay rate for even-N nuclei. Green lines represent constant separation energy 
for even-N nuclei. Faint diagonal dotted lines delineate the borders of the peak region (A — 159 to A = 167). Stable 

isotopes are shown by unfilled squares. 
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FIG. 3: Shows how the rare earth peak forms under cold environments when the r-process path encounters slower 
neutron capture rates in the peak region. Each left and right panel shows a snapshot of an abundance pattern along 
with neutron capture rates (s _1 ), /3-decay rates (s _1 ), and r-process path from a simulation with the ETFSI mass 
model. The left panel is a snapshot just after neutron exhaustion, R = 1. The right panel is a snapshot showing the 
abundance pattern and rates as the r-process path begins to move closer to stability and the peak begins to form. 
Connected black squares represent the r-process path and red lines represent constant neutron capture rates for 

even-N nuclei. All other markers are the same as Figure [2] 
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FIG. 4: A successful rare earth peak formation under a cold evolution with FRDM mass model occurs when the 
structure in the neutron capture rates lasts until the point at which /3-decays take over neutron captures in the 
region (r^r a few t^ ep ). Neutron capture rates (online red, print dotted), /3-decay rates (online blue, print 
solid), r-process path (filled triangles in the right column) and abundance snapshots (left column) are shown at 20 
(top panel), 15 (middle panel) and 10 (bottom panel) neutrons away from stability as the peak evolves during 
late-times in this cold FRDM simulation. The structure in the neutron capture rates is not yet present at 20 
neutrons away from stability (top panel). However, the structure in the neutron capture rates becomes evident at 15 
neutrons away from stability (middle panel), continuing until the abundance pattern has frozen out completely. The 

bottom panel abundances are also highlighted in Figure [l] 
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FIG. 5: A successful rare earth peak formation under a hot evolution with FRDM mass model. The r-process path 
encounters the well-defined separation energy kink in between 20 (top panel) and 15 (middle panel) neutrons away 
from stability in the peak region. The separation energy kink does not last (bottom panel) but by this time 
/3-decays have taken over and the abundances have frozen-out. Separation energies (online green, print dotted), 
/3-decay rates (online blue, print solid), r-process path (filled triangles in the right column) and abundance 
snapshots (left column) are shown at 20 (top panel), 15 (middle panel) and 10 (bottom panel) neutrons away from 

stability similar to the previous figure. 
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FIG. 6: The rare earth peak forms far from stability in this cold ETFSI simulation (top panel). However, it is 
washed out by the slow reduction in the structure of neutron capture rates of nuclei closer to stability (middle 
panel). At 10 neutrons away from stability (bottom panel) the neutron capture rates are slower than /3-decay rates, 
so that only small changes to the abundance pattern occur after this point. Hence the disappearance of the lines of 
constant neutron capture rates in the bottom panel. Neutron capture rates (online red, print dotted), /?-decay rates 
(online blue, print solid), r-process path (filled squares in the right column) and abundance snapshots (left column) 
are shown at 20 (top panel), 15 (middle panel) and 10 (bottom panel) neutrons away from stability as in previous 

figures. 
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FIG. 7: The rare earth peak begins to form far from stability in this hot ETFSI simulation (top panel) due to the 
separation energy kink structure beyond 20 neutrons from stability (not shown). However, it is washed out by the 
slow reduction in the kink structure of the separation energies of nuclei closer to stability (middle panel) . Separation 
energies (online green, print dotted), /3-decay rates (online blue, print solid), ?'-process path (filled squares in the 
right column) and abundance snapshots (left column) are shown at 20 (top panel), 15 (middle panel) and 10 
(bottom panel) neutrons away from stability as in previous figures. 
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FIG. 8: The final rare earth peak (bottom panel) is relatively flat (as compared to the solar rare earth peak) in this 

cold HFB-17 simulation. This occurs when the r-process path encounters a lack of neutron capture rate structure 
throughout the NZ-plane in the peak region. Neutron capture rates (online red, print dotted), /3-dccay rates (online 
blue, print solid) , r-process path (filled circles in the right column) and abundance snapshots (left column) are shown 
at 20 (top panel), 15 (middle panel) and 10 (bottom panel) neutrons away from stability as in previous figures. 
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FIG. 9: The final rare earth peak (bottom panel) is relatively flat (as compared to the solar rare earth peak) in this 
hot HFB-17 simulation. While there is detailed structure in the separation energies, there is little gross structure on 
the scale of the rare earth peak. The slight abundance bump off-center of the actual peak region is due to the 
complex structure in the separation energies found off-center from the peak region. Separation energies (online 
green, print dotted), /3-decay rates (online blue, print solid), r-process path (filled circles in the right column) and 
abundance snapshots (left column) are shown at 20 (top panel), 15 (middle panel) and 10 (bottom panel) neutrons 

away from stability as in previous figures. 
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FIG. 10: Highlights the nuclei which are important to rare earth peak formation. The most influential nuclei are 
closer to stability in darker shading as they set or potentially dissolve the peak structure. Also shown is the current 
extent of experimental data for the rare earth elements. Isotopes with measured masses in the AME2003 mass table 
are highlighted (online and print shaded in light gray). Recent ISOLTRAP and JFYLTRAP mass measurements are 
highlighted (online green, print dark gray) and cross section data from the online CSISRS database are highlighted 
(online red, print black). Stable isotopes are shown by unfilled squares as in previous figures. 



